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Abstract 

We present the elements of the IR-improved DGLAP-CS theory as it relates to the new MC friendly 
exponentiated scheme for precision calculation of higher order corrections to LHC physics in which IR 
singularities from both QED and QCD are canceled to all orders in a and in a s simultaneously in the 
presence of rigorous shower/ME matching. We present the first MC data comparing the implied new 
showers themselves with the standard ones using the HERWIG6.5 MC event generator as a test case 
at LHC energies. As expected, the IR-improved shower rc-populates lower values of the energy fraction 
z and lower values of the attendant pt for the standard HERWIG6.5 input parameters. Possible phe- 
nomcnological implications are discussed. 
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1 Introduction 



With the advent of the LHC, we enter the era of precision QCD, by which we mean predictions for 
QCD processes at the total theoretical precision tag of 1% or better. The attendant requirement for this 
theoretical precision is control of the 0(a 2 s L ni , a s aL ri2 , a 2 L" 3 ), n\ = 0,1,2, n-i = 1,2, 773 = 2 corrections 
in the presence of realistic parton showers, on an event-by-event basis - here, L is a generic big log. This 
is the objective of our approach to precision QCD theory, which for example will be needed for the 
expected 2% experimental precision [1] at the LHC for processes such as pp — > V + 771(7) + n (G) + X — ► 
W + ?7i' (7) + n(G) +X, V = W ± , Z, and i = e, fi, £' = is e , f M (e, fx) for V = W+(Z) respectively, and 
£ = v e , v^, t' = e,pb respectively for V = W~ . Here, we present the elements of our approach and its 
recent applications in MC event generator studies, which are still preliminary. 

At such a precision as we have as our goal, issues such as the role of QED are an integral part of 
the discussion and we deal with this by the simultaneous resummation of QED and QCD large infrared 
effects, QED ® QCD resummation [2] in the presence of parton showers, to be realized on an event- 
by-event basis by Monte Carlo methods. This is reviewed in the next Section. Let us note already that 
in Refs. [3] it has been shown that QED evolution enters at the ~ 0.3% level for parton distributions 
and that in Refs. [4] it has been shown that EW (large Sudakov logs, etc.) effects at LHC energies, as 
W's and Z's are almost massless on the TeV scale, can enter at the several % level - such corrections 
must treated systematically before any claim of 1% preciaion can be taken seriously. We are presenting a 
framework in which this can be done. The new amplitude-based resummation algebra then leads to a new 
scheme for calculating hard hadron-hadron scattering processes, IR-improved DGLAP-CS theory [5] for 
parton distributions, kernels, reduced cross sections with the appropriate shower/ME matching. This is 
summarized in Section III. In this latter Section, with an eye toward technical precision cross checks plus 
possible physical effects of heavy quark masses, we also deal with the issue of quark masses as collinear 
regulators [6-8] as an alternative [9] to the usual practice of setting all initial state quark masses to zero in 
calculating ISR (initial state radiation) effects in higher order QCD corrections. We also discuss in Section 
III the relationship between our resummation algebra and that of Refs. [10, 11], as again such comparisons 
will be necessary in assessing the ultimate theoretical precision tag. In Section IV, we illustrate recent 
results we have obtained for the effects of our new approach on the parton showers as they are generated 
with the HERWIG6.5 MC [12]. Extensions of such studies to PYTHIA [13] and MC@NLO [14] are in 
progress. Section V contains summary remarks. 

As a point of reference, in Ref. [15] it has been argued that the current state-of-the-art theoretical 
precision tag on single Z production at the LHC is (4.1 ± 0.3)% = (1.51 ± 0.75)%(QCD) © 3.79(PDF) © 
0.38 ± 0.26(£W)%, where the results of Refs. [14, 16-21] have been used in this precision tag determina- 
tionQ 



2 QED(g)QCD Resummation 

In refs. [2], we have extended the YFS theory to the simultaneous exponentiation of the large IR^ terms 
in QCD and the exact IR divergent terms in QED, so that for the prototypical subprocesses Q'Q — > 

Recently, the analogous estimate for single W production has been given in Ref. [22] - it is ~ 5.7%. 
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Q"'Q" + m(G) + 71(7) we arrive at the new result 

SUM IR (QCED) 



aa cxp = e 



^ m!n! 7 A A fc„ Al fc/ y ( 27r 

71,71=0 J Ji = l J1 32=1 32 J K 

iy(Pl+H~P2-<12~'Z2k jl ~Y1 k ' j 2 ) + D QCED 



y_ 

Y 

(i) 



P*m,n^l, ■ ■ ■ , A-m, ai , ■ ■ • , 7»<tW Q ' 

P2 12 

where the new YFS [23,24] residuals, defined in Ref. [2], fl m ,n{ki, ■ ■ ■ , k m ; k[, . . . , k' n ), with m hard 
gluons and n hard photons, represent the successive application of the YFS expansion first for QCD 
and subsequently for QED. The functions SUMir(QCED), £>qced arc determined from their analogs 
SUMi R (QCD),L»qcd in Ref. [25] via the substitutions 

&QCD — > &QCD + &QED = &QCED, 

&QCD — ► &QCD + &QED = &QCED i 

Qiils finis | rinls rinls /'0^ 

&QCD — * &QCD + &QED = ^QCED {*) 

everywhere in expressions for the latter functions given in Refs. [25] - see Ref. [2] for the details of this 
substitution. It can be readily established [2] that the QCD dominant corrections happen an order of 

magnitude earlier in time compared to those of QED so that the leading term /3q,o already gives us a 
good estimate of the size of the effects wc study. 

Important in any total theoretical precision is knowledge of possible systematic issues associated 
with ones methods. This entails the relationship between different approaches to the same classes of 
corrections and moves us to the relationship between our approach to QCD rcsummation and the more 
familiar approach in Refs. [10]. It has been shown in Ref. [26] that the latter approach is entirely equivalent 
to the approach in Refs. [11]. Establishing the relationship between our approach and that in Refs. [10] 
will then suffice to relate all three approaches. 

In Ref. [27] the more familiar rcsummation for soft gluons in Refs. [10] is applied to a general 2 — > n 

parton process [f] at hard scale Q, r x ) + f 2 (P2, r 2 ) -> h(P3, r 3 ) + fa{P4, r±) H h /„ +2 (p, l+ 2, r„ +2 ), 

where the Pi,r, label 4- momenta and color indices respectively, with all parton masses set to zero to get 



L 

C 

= J^Y. S LiH l / ] (c L ) {rth 



" c W 



where repeated indices are summed, is the jet function, Sli is the soft function which describes the 
exchange of soft gluons between the external lines, and Hf is the hard coefficient function. The attendant 
infrared and collincar poles are calculated to 2-loop order. To make contact with our approach, identify 
in Q'Q Q"'Q" + m(G) in © fx = Q, Q>, f 2 = Q', f 3 = Q", / 4 - Q'", {/ 5 , • • • , f n+2 } = {G u ■ ■ ■ , G m } 
so that n = m + 2 here. Observe the following: 



• By its definition in cq.(2.23) of Ref. [27], the anomalous dimension of the matrix Sli does not 
contain any of the diagonal effects described by our infrared functions S/jj(QCD) and Dqcd- 

• By its definition in eqs.(2.5) and (2.7) of Ref. [27], the jet function Jwl contains the exponential of 
the virtual infrared function a s $lBQCD, so that we have to take care that we do not double count 
when we use (J3)) in {1} and the equations that lead thereto. 
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It follows that, referring to our analysis in Ref. [28], we identify p*- 1 ™) in eq.(73) in this latter reference in 
our theory as 



P im) (pi,qi,P2,q2,k lr -- ,k m ) = E , . \M.y 

* — * color s. spin L 



'[/] |2 



c c , (4) 

- E *W W E E S[ f }HP(c L ) {ri} 

spnis,{ri},{r^} L — 1 L' — 1 

where here we defined J [f] = e-°- M « CD J [/1 , and we introduced the color-spin density matrix for the initial 
state, t) cs . Here, we recall (see Refs. [5,28], for example) that in our theory, we have 

e 2 as ReB QCD r ft d 3 km ft 

da n = : / I I — — ™ 6 (pi + qi - P2 - q-i - > ki) 

_(n)/ , , ,d i p 2 d A q 2 , . 

for n-gluon emission. It follows that we can repeat thus our usual steps (see Ref. [5,28]) to get the QCD 
corrections in our formula fl}, without any double counting of effects. This use of the results in Ref. [27] 
is in progress. 



3 IR-Improved DGLAP-CS Theory: Applications 

In Refs. [5,28] it has been shown that application of the result {T]) to all aspects of the standard formula 
for hard hadron-hadron scattering processes, 

o = E J dxidx 2 F t (x l )F J (x 2 )cr(x 1 X2s) (6) 

where we the {Fi(x)} and a denote the parton densities and reduced cross section respectively, leads 
one to its application to the DGLAP-CS theory itself for the kernels which govern the evolution of the 
parton densities in addition to the the implied application to the respective hard scattering reduced cross 
section. The result is a new set of IR- improved kernels [5] , 

* 2 (l-*) 7 '-/<,(7,Wl-*) , 



P qq (z) = C F F YFS ( lq )e^ 



1 - z 



in the standard notation, where 



(T) 



P Gq (z) = C F F YFS ( lq )e^ l + {l z Z) z<«, (8) 

P GG (z) = 2C G F YFS ( lG )e^{ l —z^ + T^(l ~ ^) 7G 

z 1 — z 

+ \(z 1+1G (l -z) + z(l - z) 1+ ^) - / G ( 7G )<5(1 - z)}, (9) 
P qG (z) = F YFS {lG)e^ G \{z 2 (l - z)^ + (1 - z) 2 z^}. (10) 



a s AC F 

J q = C F — (11) 



a s 4C G 

l G = C G —t=—— (13) 
So=^ + ^C4-h (14) 
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and 



FYFsilq) 



e -C E ~tq 



(15) 



r(i + 7,) 



so that 



fq{lq) 



2 2 1 



(16) 



fade) 




(17) 



C G (1 + 7g)(2 + 7g)(3 + 7g) 7g(1 + 7g)(2 + 7g) 



1 1 



(18) 



(1 + 7 g)(2 + 7g) 2(3 + 7 g)(4 + 7g) 



1 



(19) 



(2 + 7g )(3 + 7 G )(4 + 7g )' 



Here, Ce = 0.5772... is Euler's constant and T(w) is the Euler Gamma function. We see that the kernels 
are integrable at the IR endpoints and this admits a more friendly MC implementation, which is in 
progress. 

Some observations are in order. First, We note that the connection of (fT0|) with the higher order kernel 
results in Refs. [29] is immediate and has been shown in Refs. [5,28]. Second, there is no contradiction 
with the standard Wilson expansion, as the terms we resum are not in that expansion by its usual 
definition. Third, we do not change the predicted cross section: we have a new scheme such that the 
cross section in ([6]) becomes 



order by order in perturbation theory, where {P ea: P} factorize <7 U nfactorizcd to yield &' and its attendant 
parton densities {F'i}- Fourth, when one solves for the effects of the exponentiation in (fT0|) on the actual 
evolution of the parton densities from the typical reference scale of Qo ~ 2GcV to Q = 100 GeV one 
finds [5,28] shifts of ~ 5% for the NS n=2 moment for example, which is thus of some phenomcnological 
interest- see for example Ref. [30]. Finally, we note that we have used [2] the result ^ for single Z 
production with leptonic decay at the LHC (and at FNAL) to focus on the ISR alone, for definiteness 
and we find agreement with the literature in Refs. [31-33] for exact 0(a) results and Refs. [34-36] for 
exact 0{a 2 s ) results, with a threshold QED effect of 0.3%, similar to that found for the parton evolution 
itself from QED in Refs. [3]. Evidently, any 1% percision tag must account for all such effects. 

3.1 Shower/ME Matching 

In using {J) in (|20p for a'faiXj), we intend to combine our exact extended YFS calculus with HERWIG [12] 
and PYTHIA [13] as follows: they generate a parton shower starting from (x\,X2) at the factorization 
scale fx after this point is provided by the {F[} and we may use [2] either a pT-matching scheme or a 

shower-subtracted residual scheme where the respective new residuals {/3„ jm (fci, . . . , fc n ; k[, . . . , k' m )} are 
obtained by expanding the shower formula and the result in |T|) on product and requiring the agreement 
with exact results to the specified order@ This combination of theoretical constructs can be systematically 
improved with exact results order-by-order in a s ,a, with exact phase space. |f| The recently developed 
new parton evolution algorithms in Refs. [38] may also be used here. 

The issue of the non-zero quark masses in the initial state radiation is present when one wants 1% 
precision, as we know that the parton densities for the heavy quarks are all different and the generic size 

2 See Ref. [37] for a realization of the shower subtracted residual scheme in the context of QED parton showers. 
3 The current state of the art for such shower/ME matching is given in Refs. [14], which realizes exactness at 0(a s ). 




(20) 
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Histogram of EF of parton shower constituents in herwig6.5 for QCD 2->2 hard parton scattering. 
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Figure 1: The z-distribution shower comparison in HERWIG6.5 - preliminary results. 



Histogram of EF of parton shower constituents in herwig6.5 for QCD 2->2 hard parton scattering. 



DGLAP-CS 
IR-lmp.DGLAP-CS 



100000 - 



80000 - 




Figure 2: The z-distribution shower comparison in HERWIG6.5 at small z - preliminary results. 



of mass corrections for bremsstrahlung is a s /7r for cross sections [39], so that one would like to know 
whether regularizing a zero-mass ISR radiation result with dimensional methods, carrying through the 
factorization procedure gives the same result as doing the same calculation with the physical, non-zero 
mass of the quark and again carrying through the factorization procedure to the accuracy a 2 /ir 2 , for 
example. Until the analysis in Rcf. [9], this cross check was not possible because in Refs. [6,7] it was 
shown that there is a lack of Bloch-Nordsieck cancellation in the ISR at 0(a 2 ) unless the radiating quarks 
are masslcss. The QCD rcsummation algebra, as used in (QJ, allows us to obviate [9] this theorem, so 
that now such cross checks are possible and they are in progress. 



3.2 Sample MC data: IR-Improved Kernels in HERWIG6.5 

We have preliminary results on IR- improved showers in HERWIG6.5: we compare the z-distributions 
and the pt of the IR-improved and usual DGLAP-CS showers in the Figs. 1, 2, 3. As we would expect, 
the IR-improved shower re-populates the soft region in both variables. The details of the implementation 
procedure and the respective new version of HERWIG6.5, HERWIG6.5-YFS, will appear elsewhere [40]. 
The analogous implementations in PYTHIA and MC@NLO are in progress, as are comparisons with IR 
safe observables. 
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Histogram of PT 2 for QCD parton shower in Herwig6.5 for 2->2 hard parton scattering. 



DGLAP-CS 
IR-lmp.DGLAP-CS 




Figure 3: The p-r-distribution shower comparison in HERWIG6.5 - preliminary results. 

4 Conclusions 



The theory of Ref. [23] extends to the joint resummation of QED and QCD with proper shower/ME 
matching built-in. For the simultaneous QED®QCD resummed theory, full MC event generator realiza- 
tion is open: a firm basis for the complete O(o^,aa s ,a 2 ) MC results needed for precision LHC physics 
has been demonstrated and all the latter are in progress - see Refs. [41] for new results on e expansions 
for the higher order Feynman integrals needed to isolate the residuals in our approach for example. This 
allows cross check between residuals isolated with the quark masses as regulators, something now allowed 
by the result in Ref. [9], and those isolated in dimensional regularization for the massless quark limit. 
Such cross checks are relevant for precision QCD theory. The first MC data have been shown with IR- 
improved showers in HERWIG6.5. The spectra are softer as expected. We look forward to the detailed 
comparison with IR safe observables as generated with IR-improved and with the usual showers - this 
will appear elsewhere. [40]. Already, semi-analytical results at the~/3o'o arc consistent with the literature 
on single Z production, while a cross check for the analogous W production is near. As the QED is at 
0.3% at threshold, it is needed for 1% precision. 
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